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We present spectral and optical properties of the Hubbard model on a two-dimensional square 
lattice using a generalization of dynamical mean-field theory to magnetic states in finite dimension. 
The self-energy includes the effect of spin fluctuations and screening of the Coulomb interaction 
due to particle-particle scattering. At half-filling the quasiparticles reduce the width of the Mott- 
Ifubbard 'gap' and have dispersions and spectral weights that agree remarkably well with quantum 
Monte Carlo and exact diagonalization calculations. Away from half-filling we consider incommen- 
surate magnetic order with a varying local spin direction, and derive the photoemission and optical 
spectra. The incommensurate magnetic order leads to a pseudogap which opens at the Fermi energy 
and coexists with a large Mott-Hubbard gap. The quasiparticle states survive in the doped systems, 
but their dispersion is modified with the doping and a rigid band picture does not apply. Spectral 
weight in the optical conductivity is transferred to lower energies and the Drude weight increases 
linearly with increasing doping. We show that incommensurate magnetic order leads also to mid-gap 
states in the optical spectra and to decreased scattering rates in the transport processes, in quali- 
tative agreement with the experimental observations in doped systems. The gradual disappearence 
of the spiral magnetic order and the vanishing pseudogap with increasing temperature is found to 
be responsible for the linear resistivity. We discuss the possible reasons why these results may only 
partially explain the features observed in the optical spectra of high temperature superconductors. 

PACS numbers: 71.27.-)-a, 74.72.-h, 75.10.-b, 79.60.-i. 



I. INTRODUCTION 

In the past decade the interest in the physical prop- 
erties of correlated electronic systems has greatly in- 
creased. One reason is the strong local correlations on 
transition metal ions in cuprate superconductors and 
manganites, and the corresponding unusual properties of 
these compounds. The parent undoped compounds are 
Mott-Hubbard or charge transfer insulators, while dop- 
ing leads to correlated metals in which the kinetic cru 
ergy of charge carriers competes with magnetic order.EJ 
One of the most spectacular consequences is the onset 
of high-temperature superconductivity in the cuprates. 
It is believed that a satisfactory description of the nor- 
mal phase properties is a prerequisite for the understand- 
ing of the microscopic mechanism of pairing in high- 
temperature superconductors. The electronic states in 
Cu02 planes of cuprate superconductors are usually de- 
scribed in terms of the Emery model which includes 
the hvbridization between Cu(3da;2_y2) and 0(2px(y)) 
states.D However, hole doging leads to the formation of 
local Zhang-Rice singletsja and the essential excitations 
in the cuprates within a window of a few eV around the 



chemical potential are well reproduced using the effec- 
tive two-dimensional (2D) Hubbard model with extended 
hoppingD and a large local Coulomb interaction U , as 
shown by various numerical studies of the t-J and Hub- 
bard modelE 



Recently it was showroiJ'lil that the effective Hub- 
bard model has to include hopping beyond the nearest- 
neighbors. The second nearest-neighbor hopping changes 
the dispersion of the quasiparticle (QP) states, and is 
therefore crucial for understanding angular resolved pho- 
toemission (ARPES) data of the antiferromagnetic (AF) 
insulator Sr2Cu02Cl2.u Both second and third neighbor 
hopping parameters follow from the down-folding proce- 
dure in electronic structure calculations |3 and influence 
the shape of the Fermi surface. They have a particular 
relation to the value of the sunerconducting transition 
temperature at optimal doping.Q 

The superconductivity occurs in the cuprates under 
doping 5 = 1 — n of a half-filled (n = 1) AF insula- 
tor, and is accompanied by a gradual modification of 
the magnetic order. The nature of magnetic correla- 
tions in doped materials is therefore a central issue in 
the theory of the cuprate superconductors. Undoped 
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La2Cu04 is a commensurate AF insulator, while doping 
by Sr into La2-a;Sr2:Cu04 results in short-ranaep4F or- 
der within incommensurate magnetic structures£H13 Such 
an incomm en| S ur ate magnetic order waa-.iiwjiced found 
analyticallyP'ta in HartreeyEdck (HF),ytj and in a 
slave-boson approximation.Ej~t3 However, in order to un- 
derstand the transport properties one has to go beyond 
an effective single-particle description and include the dy- 
namics due to local electron correlations. 

A sufficiently accurate treatment of local electron cor- 
relations remains one of the challenging problems in mod- 
ern solid state theory. Although an important progress in 
the present understanding of strongly correlated fermion 
systems occurred recently due to numerical methods, 
such as quantum Monte Carlo (QMC) and exact diago- 
nalization (ED), an analytic treatment that maintains lo- 
cal correlations is needed to investigate the consequences 
of strong correlations in the thermodynamic limit. An at- 
tractive possibility is the limit of large spatial dimension 
(rf — cxd), where the diagrams in the perturbative expan- 
sion collapse to a single site and-the fermion dynamics is 
described by a local self-energy r3 This allows a mapping 
of lattice models onto quantum impurity models, which 
can then be solved self-consiatently using the dynamical 
mean-field theory (DMFT).^ 

The DMFT was quite successful for the Hubbard 
model with nearest-neighbor hopping t at half-filling, 
where it predicts the Mott transition to the insulating 
state (n = l).c2l|TIhis was also found by Logan, East- 
wood, and TuscliEiJ for the d = oo case using an ana- 
lytic method. Attempts to use DMFT at arbitrary filling, 
however, made it clear that the local self-energy becomes 
particularly important in systems with magnetic long- 
range order (LRO), which are easily destabilized when 
the correlation effects are overestimated. The self-energy 
therefore plays a decisive role and has to be described. 
beyond second order perturbation theory (S0PT).E3'E^ 
This has made the application of DMFT to magnetically 
ordered systems notoriously difhcult. Recently we have 
shown that the screening of local Coulomb interaction 
by the particle-particle diagrams plays a crucial role in 
stabiliziiig the incommensurate magnetic LRO in doped 
systems 

The advantage of using DMFT becomes clear by look- 
ing at the single hole problem, which can be solved ex- 
actly in the d —> oo limit .E3 The method becomes exact 
because the quantum fluctuations are of higher order in 
the 1 /d-expansion than the leading potential term which 
originates from the Ising part of the superexchange in- 
teraction J = At'^/U. Therefore, applying DMFT to the 
d = 2 case might still capture the essential features that 
result from the coupling of a moving hole to local spin 
fluctuations. We will show below that in fact such quan- 
tities as the spectral function, the QP band, and the size 
of the QP spectral weight are well reproduced within the 
DMFT, which for a single hole includes only those pro- 
cesses which are present in the t — model. Although 
this approach becomes exact only in the d ^ oo limit ,ca 



it gives a sufficient accuracy of the onfeBarticle spectral 
function even in finite dimension d = 2.Ej The DMFT al- 
lows us to calculate the optical ccmductivity in the d = oo 
limit of Metzner and Vollhardtt3 from the knowledge pi 
the local self-energy without further approximations.EZI 
The studies performed in this limit for the nonmagnetic 
systems already allowed a qualitative reproduction of 
such experimental observations in the cuprates as the in- 
crease of the Drude peak with doping, ancLaiemperature 
and doping dependent mid-infrared peak.oEj 

The paper is organized as follows. The self-consistent 
procedure to determine a local self-energy within the 
DMFT is introduced in Sec. II. It consists of the HF 
potential and the dynamical part due to spin fluctua- 
tions which uses a Coulomb interaction rcnormalized by 
particle-particle scattering. The formalism to calculate 
the one-particle and optical excitation spectra in the spin 
spiral (SS) states is developed in Sec. III. Next we an- 
alyze the numerical results for the one-particle spectral 
properties in Sec. IV, where we show how they change 
with doping and with increasing temperature. The opti- 
cal properties are presented in Sec. V; there we discuss 
the new effects in the optical conductivity, scattering rate 
and the effective mass which arise due to extended hop- 
ping and by increasing the value of Coulomb interaction 
U. Sec. VI presents a short summary and conclusions. 



II. DYNAMICAL MEAN-FIELD THEORY FOR 
SPIN SPIRAL ORDER 

A. Dynamical mean-field equations 

We consider the spectral and optical properties of the 
minimal model for strongly correlated electrons in high- 
temperature supeEcanductors, the Hubbard model with 
extended hopping, at] 

H = -Y^t,jal^aj^ + UY^n,fn,i, (2.1) 

ija i 

where aj^ is a creation operator of an electron with spin 

a at site i, and riio- = a\^a^^. The hopping elements 
tij — t, t' and t" stand for the nearest neighbor, second- 
nearest neighbor, and third-neighbor hopping on a 2D 
square lattice, and serve to model the electronic states of 
the charge-transfer type in the cuprates. For convenience 
we choose t as the energy unit. 

It is interesting to note that hopping beyond nearest- 
neighbors contributes to the energy and other properties 
not only in a 2D model, but also in the limit of c? — s- oo. 
The energy contributions due to more distant neighbors 
are finite due to the scaling of the hopping parameters 
on a hypercubic lattice. It is given by ~ (i-ll*-ill/2 
(see Refs. p^pO| ), where ||i — j|| is the distance between 
i and j defined by the 'bond metric', and gives the scal- 
ing factors ^ 1/Vd for first-, and ^ l/d for second- and 
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third-nearest neighbor hopping, respectively, as the latter 
sites are two bonds apart. 

As mentioned above, we adopt the limit of infinite di- 
mensions to determine the spectral properties of the Hub- 
bard Hamiltonian on a square lattice in the thermody- 
namic limit. In order to simplify the numerical evaluation 
of the self-energy we introduce an ansatz for the modified 
magnetic order in the doped systems, and assume incom- 
mensurate SS structures with a large but finite period- 
icity. This approach captures the essence of the compe- 
tition between the weakened short-range AE order, and 
the kinetic energy induced by hole doping,ll3 and allows 
to treat the systems in the thermodynamic limit at low 
temperature. The spiral states are characterized by the 
amplitude of the local magnetization. 



mo = - nil) I, 



(2.2) 



which is independent of the site index i. The direc- 
tion of the magnetic moment at each site i is specified 
in the global reference frame by two spherical angles, 
= (0i:^i)i a-iid therefore the original fermion oper- 
ators, {aj|,aj|}, are transformed to the new fermions 
quantized-^ith respect to the local quantization axis at 
each site,li3 



Act ' 



(2.3) 



where 7^(^^^) = e-'('^'/2)*=e-'(^'/2)*f is the rotation ma- 
trix, (jji and 9i are polar and azimuthal angle, respec- 
tively, and (7y and are Pauli spin matrices. This trans- 
forms the Hubbard Hamiltonian ( |2.1[ ) to the following 
form, 



H=- ^^,,cL[7^t(^^,)7^(^■; 



(2.4) 



In the SS states we take the polar angle site-independent, 
Oi = 9, and the azimuthal angle is given by the wave- 
vector Q of the spiral, 0i = Q • Ri- Using the period- 
icity of the Ti^TZ matrix in Eq. {2A), one finds after a 



Fourier transformation that the kinetic energy takes a 
simple (2 x 2) matrix form, 

+ 5 £k+.s(i - COS0 (7^ + Sin0 (7a;), (2.5) 

where Ck = —2t{coskx -f cosfcj,) — 4i' cos /c^: cos fcj^ — 
2<"(cos2fca; +cos2ky) is the electron dispersion in a non- 
interacting system in the global reference frame. Here 
we limit ourselves to plane spirals, and choose 9 = 7r/2. 
Therefore, the order parameter rotates in the (a, h) plane, 
(S,) = (mo/2)[cos(QRj),sin(QR,),0]. Double spirals 
were shown to be unstable in the 2D t-J model,E3 and we 
have no reason to believe that they might be stabilized 
by further neighbor hopping. 



In order to construct the leading local part of the self- 
energy, we use the DMFT and consider the impurity 
model coupled to the lattice by the effective field (for 
more details see Ref. ^9|). The Anderson model of a 
magnetic impurity coupled to a conduction band with 
SS order consists of a "non-degenerate impurity orbital" 
at site o, with the fermion operators {ftm focr\^ ^'^d the 
conduction electron bath as an "effective SS conduction 
band" described by the operators {c\^„,c^cr'}^ 



H.mp^ Sf flf,, + [TQ(k) 



+ E [foa [^Q(k)]CTCT'CkCT' + H.C.] + Uni^ni^, (2.6) 



kCTCT' 



where is an impurity energy level, and TQ(k) is an 
effective o ne-p article energy of the same functional form 
as rQ(k) { p,.5\ j. The hybridization (2 x 2) matrix in the 
local reference frame. 



1 

2' 



cos-(l + a^) - sin -(7+ 

9,- , 9 
cos -(1 - <7z)+ 



(2.7) 



where a± = i id'y^ is given by the individual hy- 
bridization elements in the glob al r eference frame, Wk = 
e^'^'^'Uoi. The Hamiltonian ( ^.6| ) is quadratic in Cj.^^, 
and the bath of conduction electrons can be integrated 
out giving raise to an effective actkui of the impurity 
electrons which is of the usual form,ta 

„. Jo 



+ U I drnf^(r)nf^(T), 



(2.8) 



where are Grassmann variables for the /- 

electrons. The Weiss effective field GQa-a'i''' ~ ''"') ^ 
(2 X 2) matrix in spin space, 



E 

k 



Fp(k)[*c..-rQ(k)]-if^(k). (2.9) 



For a plane spin spiral with 9 = ti/2 the Weiss effective 
field becomes a diagonal matrix in spin space. 



^QCTCT'(*^n) ~ So 



(2.10) 



Note that this resu lt only depends on the functional form 
of Eqs. (2.5)-(2.7), and not on the parameters, except 
that it holds for a plane spiral. It implies that the lo- 
cal spin fluctuations are decoupled from the local charge 
fluctuations, and simplifies the present self-consistent cal- 
culation for SS states within the DMFT approach, as all 
local quantities including the self-energy E, are diagonal. 
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In the spirit of the DMFT approach, we acproximate 
the Green function using a local self-energy 



GQ^(k, iLUi,) 



Q i^^i^), 

(2.11) 



where /i = — £/ is the chemical potential. The lattice (fi- 
nite) dimensionality enters via the one-particle energies 
TQ(k) and gives rise to the k-dependence of the spec- 
t ral f unction. The lattice one-particle Green function 
( ^.11 ) is described by a (2 x 2) matrix G'Q(k, iw,y) in 
spin space, where Ui, are fermionic Matsubara frequen- 
cies. The corresponding local lattice Green function, 
GQ(iWy) = iV"^ GQ(k, zwi.) cx Jo-o-', is diagonal in 
spin space due to the parity of the kinetic energy Tq (k) . 

The self-energy consists of the HF part, Sq^ = U{nog-) 
with a — — (T, and the spin-fluctuation (SF) part, 
I]Q^(icj^), which is determined by the many-body ef- 
fects. Using the cavity met/iod@ for a hypercubic lattice 
at d = oo, we verified that the dynamical Weiss field, 
C/q^(za;y), can be computed from the Dyson equation 
of the Anderson impurity model (2.6) with broken spin- 
symmetry, 

g^aituj,) -1 - GQituj,,) -1 + t^f^ituj,) . (2.12) 

Eqs. (|2.1l|) and ( |2.12|) are fundamental in the DMFT,[i 
and can be solved self-consistently, provided an expres- 
sion for the self-energy is known. 



B. Thermodynamic potential at finite temperature 

The calculations at finite temperature T require the 
knowledge of the free energy, F{T, N^), being a thermo- 
dynamic potential for a system of = Nn electrons. It 
has to be minimized to find a stable SS state which de- 
termines the system properties. The free energy may be 
found from the grand canonical potential, ^(T, fi), using 
the standard approach for quantum many-body systems, 
F{T,Ne) = ^iT,iJ,) + /xiVeB For a translationally in- 
variant lattice model with local self-energy one finds the 
functional form of the grand canonical potential, 

i](T,A*) = r!o(r,/i) + <i>sF[GQ] 

-/3-i^lndet \l-G'l^{k,iuj^)±^{ioj^) 



(2.13) 



with fi = l/fcsT, and Sq (icji,) is the self-energy dis- 



cussed below. The functional ( 2.13| ) is stationary, i.e., 
5VL = ensures that the minimum of the grand canonical 
potential has been found, and determines the self-energy 
from the Luttinger-Ward functional. 



G, 



QJ 



N 



(2.14) 

Our perturbative expansion is constructed around the 
symmetry-broken HF state, hence the grand canonical 
potential of the "non-interacting" reference system in- 
cludes a correction term to avoid double-counting and 
reads. 



^o{T,^J) = /3-i^lndet 



G°j(k, lu. 



-UN{no^){noi). (2.15) 

The spectrum which defines flo{T,fi) is given by the 
Green's function in the HF approximation, 

-1 _ „■, 



G2j(k, iLo,) 



(2.16) 



The Luttinger-Ward functional $^^[Gq] in Eq. (|2.13| ) 
is defined via the diagrammatic expansion of Sq m 
terms of the full Green's function Gq. The self-energy of 
the infinite-dimensional Hubbard model is a local dynam- 
ical quantity and involves only the local component of 
the Green 'ii_function ( 2.11 ). This imphes that <i>^F[G] = 
^*fmp[G]J^ meaning that the functional $^^[G] can be 
approximated by some infinite subset of the one-particle 
irreducible closed Feynman diagrams of the Anderson im- 
purity model (2.6). We-|take for ^^^[G] the sum of all 
particle- hole diagrams |2j and the effective particle-hole 
interaction C/,oE3 

SP-iV(02 + 0± + 0||), (2.17) 
l._,,.o^ fa,,. . mu. . ^2.18) 



0± 



In 1 



E (*^m)Xx7 (it^p), 



(2.19) 



2 



(2.20) 



where 



(2.21) 

is the noninteracting particle-hole susceptibility. Self- 
consistency would require that = $^F[Gq]; here 
instead we apply the non-self-consistent procedure intro- 
duced by Bulut, Scalapino, and White,Ej and approxi- 
mate $^^[Gq] $^^[^q]- It has been shown that this 
procedure may be regarded to be a reas onab le approxi- 
mation as the thermodynamic potential ( ^.13 ) is station- 
ary and one expects not to move too far away from its 
minimum. 
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C. Self-energy with local spin fluctuation 



U = U/[l + Ux''PiO)i 



(2.25) 



It is known to be notoriously difficult to obtain an an- 
alytic expression for the self-energy, and so far mostly 
an ansatz within the iterative ppfiturbation scheme (IPS) 
based on SOPT has been usedE3 The ansatz introduces 
an approximate form of self-energy which starts from the 
SOPT and allows to reproduce the exact results in cer- 
tain limits. Although this approach reproduces the cor- 
rect large U limit ,E3 it overestimates the correlation ef- 
fects in the nonmagnetic states, and thus becomes un- 
controllable in the intermediate U regime. Therefore, it 
cannot be applied to investigate the phase stability and 
dynamics in the magnetic states of the Hubbard model. 
We have verified that the AF LRO disappears in the 2D 
Hubbard model {f = t" = 0) at half-filling for U ~ 5t for 
t' = t"j-^ 0, if the formul a intr oduced by Kajueter and 
KothaiEa is used (see Sec. |lVA|) . 

The SF part of the self-energy, Il^^{iuj^) follows from 
the Kadanoff-Baym potential ( 2.13| ) containing a class of 
diagrams up to infinite order, 

Here we approximated T,[G] by and avoid self- 

consistency. The transverse part in Eq. ( 2.22 ) resembles 
the self-energy derived by the coupling of the moving 
hole to transverse-spin-fluctuations, as derived using the 
spin- wave theorytZI However, the longitudinal part is not 
included in the latter approach, and we find that it can- 
not be neglected in the relevant regime of parameters for 
high temperature superconductors. 

The self-energy in a magnet ic sy stem is calculated 
using the Weiss effective field (2.12|) in the symmetry- 



broken magnetic state. The transverse. 



UxTa^^,) 



and longitudinal. 



^(0) 



1 - C/2xj°^(icJp)X|°^(iwA') 



(2.23) 



(2.24) 



susceptibility in Eq. ( |2.22 ) are found in random phase 
approximation (RPA) with renormalized interaction U . 
Here the non- interacting susceptibilities, ^^'^(ici;^), are 
calcu lated from the dynamical Weiss field Green func- 
tion ( |2.12| ). 

We would like to emphasize that the, renormalized in- 
teraction tj is not a fitting parameter J2j but results from 
the statieirSfireening by particle-particle diagrams which 
leads tocaEa 



where the particle-particle vertex is again determined by 
the Weiss field, 

X^^(O) -r'E^QT(*^M)^QiH^A.)- (2-26) 



This pScreening effect giv-es the magnetic str ucture 
factoiO and the self-energyEj calculated from Eq. ( 2.22| ) 
in good agreement with the QMC results, and depends 
on the underlying magnetic order. It is largest in the 
paramagnetic states and vanishes in the Neel state at 
n = 1 for [/ — !■ cx), and is thus very imppctant when the 
magnetic p hase d iagrams are considered.E3 The proposed 
self-energyi 2.22| ) expresses the spin- fluctuation exchange 
intcractioncS with m effective potential due to particle- 



particle scattering.ca 

Eqs. (2.11), ( ^.12 ) and ( 2.22| ) represent a solution for 
the one-particle Green function within the DMFT. They 
have been solved self-consistently and the energetically 
stable spiral configuration was fowid. This procedure is 
further justifled by the sum rule,Ej 



^^E,(zc..)G.(z...)e-"°' 



[/(riotnoi), (2.27) 



which is well fulfilled in --the present approach with 
Ujnq ^noi) ~ C7(no|) (rioi)-E^ W e also show below (Sec. 
IV A| ) that the local self-energy ( ^.22 ) leads to an overall 
satisfactory agreement with the QMC and ED data. 



III. EXCITATION SPECTRA 
A. Photoemission at finite temperature 

A complete theory of photoemission (PES) would re- 
quire an analysis not only of the one-particle Green's 
function but also of the three-particle Green's functions. 
We would like to point out that quantitative calculations 
of the three-particle Green's function for strongly cor- 
related systems has not yet proven feasible. However, 
some aspects of the problem can be discussed in terms of 
the one-electron spectrum, provided that "flnal-state" or 
"particle-hole" interactions can be neglected. Under this 
assumption the problem simplifies and the PES spectrum 
may be determined using the one-particle Green function 
alone. Such an approach which neglects final-state and 
particle-hole, interactions has been applied with success 
to interpretcl the diS|]pr.rsion found in the ARPES data of 
the copper oxides .CjO 

Here we shall derive the relation of the PES spectra 
to the one-electron spectral function within the "sud- 
den" approximation, where final state interactions are 
neglected.c3. To be specific, let us consider a transition 
from a state \n) with energy En into a state of the form 
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^ki/l™)' which we treat the photoelectron in state 
[kv) as dynamically decoupled, but retain the full many- 
body interactions of the electrons in the bulk described 
by the Hubbard model Hamiltonian. The PES intensity 
for the magnetic system with an incommensurate mag- 
netic order is nontrivial within the DMFT approach, as 
one cannot use a Bogoliubov transformation to establish 
the relation between the measured electrons and their 
local states in the SS state. The outgoing photoelectron 
is observed in the global reference system, whereas the 
quantum states of the bulk \n) have to be considered 
within the local reference system for the spin degrees of 
freedom. For clarity we write in the following the oper- 
ators for the scattered states in capital letters and the 
operators for the electronic states of the solid described 
by the present model Hamiltonian in lower-case letters. 

At finite temperature T we consider the probability 
density of the absorption of a photon with frequency uj 
in a grand canonical ensemble and obtain, 

mn 

(3.1) 

with Kn = En — and the partition function Z — 
e~^^" . The amplitude of the transition, 

= (m|^k. ^p..[^q(P' p')]..' v^s,. I«> (3-2) 



is determined by the optical matrix element, 

'•19 9 
Aq(p,p') = -Ap,p4cos-(i - <T^) +sin-(7_] 

1. 9 9 

+ -Ap,p,_Q[cos-(i +0-^) -sin-(7+], (3.3) 

and can be calculated using the Bloch wave-functions in 
the global reference system, 



(3.4) 



where e is the polarization vector. The operator k = 
— iV conserves total momentum in the scattering plane, 
so that Ap pi (X 5p||,p'+K, where K is a 2D lattice vec- 
tor, and p|| is the photoelectron momentum component 
in the 2D plane. 

For solids the outgoing wave solution is the "tirrie=. 
inverted low-energy electron diffraction (LEED) state" £3 
The LEED state consists of an incoming plane wave, re- 
flected plane waves, and a combination of Bloch waves in- 
side the solid which fulfill the matching boundary condi- 
tions. In lowest order we have one (damped) Bloch wave 
travelling away from the surface. In the time inverted 
(complex conjugated) state the Bloch wave travels to- 
wards the surface, and goes over in a plane wave outside. 
The LEED scattered waves become incoming waves on 



time inversion, and give no contribution to the photocur- 
rent. The photoelectron is usually detected at energies 
which are much higher than the typical energy regime de- 
scribed by the Hubbard model, and therefore the Bloch 
waves occupy high-energy quantum states which are ini- 
tially unoccupied, 

A. \n) ^ 0. (3.5) 
Hence, we obtain for the plane SS state {9 — 7r/2), 

W(k, c^) = -i V I Ak.k|'nF(ek - c^) 

TT ^ — ' 

xImC-(k-Q/2,£k-c^), (3.6) 

where nF{(^) is the Fermi function, and the following 
identity, vahd only for plane spirals {9 — Tr/2), has been 
used (Ao- = 1,-1 for a =t, i). 



E 



Gaa' (k - Q/2, u;) = Yl G--' + 



(3.7) 

Within the "sudden" approximation the measured 
PES spectra near the Fermi energy can therefore be re- 
lated to the one-electron spectral function Eq. (3.6) of 
the system with local spin-quantization axes, defined by 



ImGo-CT' (k, w) = - 



^(n|4^,|TO)(m|ck^|n) 



xe-'^'^-Siiv- Kn + K,n). (3.. 



Therefore, the total one-particle excitation spectra is de- 
scribed by the spectral function. 



A{k,uj) = - -V ImG, 

TT ' ^ 



Q.trtr' I k - + ie ) , (3.9) 



where GQ^aa' (k— Q/2, w-l-ie) is given by Eq. ( ^.11 ), and 
a numerical broadening e > 0. The electron occupation 
number (rik) normalized per one spin, equal to the one- 
electron removal sum, can be obtained without analytic 
continuation of the Matsubara Green's function (2.11) 
by performing a direct summation over the Matsubara 
frequencies. 



^'^k) = ^ E e^"''°'GQ..,, (k . (3.10) 

Finally, we calculate also the total densities of states in 
the AF and SS states using the derived spectral functions 



(3.11) 
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B. Optical conductivity 



kx{i - I) 



El 



IWlt, 



We derived the complex optical conductivity (T^xi^) 
for the spiral magnetic order following the formalism in- 
troduced by Shastpjinand Sutherlandjlj and by Scalapino, 
White and Zhang.E3 Their derivation has to be general- 
ized to the case of extended hopping. Moreover, as the 
symmetry is locally broken in a magnetic system with lo- 
cal quantization axes, the calculation of the optical con- 
ductivity is not straightforward. The Hubbard Hamilto- 
nian (2.4) within the local reference system for the spin 
quantization axis and first-, second-, and third-nearest 
neighbor hopping elements, tn = t,t' , t" , respectively, has 
an electron kinetic energy, 



il.(7(7' 



(3.12) 



where ' indicates a restricted sum, with R; = + 
\\i — /||a;X -|- \\i — l\\yy around a given lattice site i, and 
X — (1,0), y — (0,1) are unit lattice vectors. We in- 
troduce a directed "bond metric", \\i — l\\x{y), which is a 
distance between two sites, i and I, on the lattice, and 
counts the number of a;(j/)— oriented bonds that connect 
site i with site I, respectively, e.g., \\i — l\\x = 2 and 
||i — Z||y = if the electron hops to a third- nearest neigh- 
bor with amplitude t" along an a;— oriented link. Here 
TZ{fli) is the unitary matrix which transforms the origi- 
nal fermions, {al^ '^Ixi' ^^^^ fermions quantized with 
respect to local quantization axes at each site, {c|p c||}. 



introduced in Eq. (2.3). In what follows we are interested 
in the current response to a vector potential along the x- 
direction of the 2D square lattice Ax{l, t). In the presence 
of a vector potentialpthe hopping term is modified by the 
Peierls phase factor ,c£l either e-Kp{+ieAx{l,t)\\i — l\\x} or 
e'icp{—ieAxil,t)\\i — l\\x}^ for tu or tu, respectively. Ex- 
panding these phase factors in the usual manner up to 
second order ~ one finds. 



e£{i-l)Ax{l) + -kx{i~l)Ax{lf 



(3.13) 



Here jxii ~ the x-component of the paramagnetic 
current density. 



3 



-i) = iJ2\\^- iWxtu \cl [7^t(^7,)7^(^7^)]^^, 



- cl, [7^t(^^07^(a)].vC-] ' (3-i4) 

and kx(i — l) is the kinetic-energy contribution due to the 
a;— oriented links, weighted by the metric factor connect- 
ing site i with site 



+ cl,[n*{n,)n{n,)l,,c,,\, (3.15) 

After performing a Fourier transformat ion o ne finds 
the average contribution of kinetic energy (3.15) per one 
site. 



•kcr'/i 



(3.16) 



k.crcr' 



with the coupling between the transformed elements at 
momenta k— Q/2 and k+Q/2 due to the magnetic order. 



4,Q(k) = ^ex(k~ -y ) (i -l-o-^cosf 



ax sm f 



\ex(k + - a^^cose + ax sin 6*) , (3.17) 



and ea;(k) = —2t cos kx — 4:t' cos kx cos ky — 8t" cos 2kx. 

As usually, the optical conductivity in long-wavelength 
limit q — > 0, axx{i^) = o'xxi^) + i'^xxi^)^ determined 
by the cup^nt response to a vector potential along the a;- 
direction,E3 and one finds using the Kubo linear response 
theory. 



CTxxi^) = -e 



2 \- kx) - J^xxiq = 0,LU + iO+) 

i{uj + iO+) 



(3.18) 



where Axxi^., i^^j) is the current-current correlation fimc- 
tion. 



1 

Axx{ci, iujf,) = — J dTe''^>^''(jxici,T)jx{~q, 



0)). (3.19) 



The latter correlation functiaO|_is given exactly by the 
particle-hole bubble diagram£Zlo where for q ^ 0, 



(3.20) 



and for the present SS state. 

If Q\ - 
jx,Q(k) ^ - jxlk - — j {1 + azCos9 - ax sin 9) 

+ ^ 3x (k + § ) (i - cos 9 + ax sin 9), (3.21) 

with jx (k) = 2t sin kx + 4t' sin kx cos ky + At" sin 2kx ■ The 
advantage of using the DMFT with the local self-energy is 
that the vertex corrections to the current-current correla- 



tion function ( 3. IE ) disappear, and the optical conductw- 
ity can be calculated without further approximations.E3 
We have verified for large variety of doping levels and 
temperatures that the following optical sum rule. 



2 / dcja'^^iuj) = e^TT {-kx), 
lo 



(3.22) 
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is always fulfilled within the numerical accuracy, in con- 
trast to the approaches whic h cannot be derived in a 
diagrammatic way. Eq. ( 3.22| ) is also used to define the 
plasma- frequency ujp, 



= 8 / dcjcr^^(w). 



(3.23) 



For the discussion of the complex conductivity function, 
it is convenient to introduce the following parametriza- 
tion by the scattering rate, T(tj)^^ , and,the effective mass 
m*{uj)/me {irie is the electron mass),E3 



47r 



1 



^[Lj) — ILJ ^ 



(3.24) 



^From the real part of the optical conductivity ( 3.18| ) 
we find in the limit w — > the static conductivity, 

(T^^(w = 0) = e^TrD -I- lim -ImA^^iq = 0, w), 

(3.25) 

with the Drude weight D which may be obtained from the 
zero-tcmpcrature extrapolation of the currcntpCjarrent 
correlation function in the upper complex planejl3 



D 



lim [ 



Re A^:,(q = 0,27riT)]. (3.26) 



The optical conductivity allows to determine the in-plane 
static resistivity. 



(3.27) 



where the static conductivity cr'^x{uj = 0, T)~^ is obtained 
as in Eq. ( ^.25| ). We present the results obtained for the 
optical conductivity and static resistivity in Sec. 0, and 
show that the magnetic order in the doped compounds 
has directly measurable consequences for these quanti- 
ties. 



IV. ONE-PARTICLE SPECTRA 



order is there of crucial importance. Therefore, the at- 
tempts to describe the AF order based on the SOFT 
within the lES failed and the magnetic order disappeared 
at larger [/.l3 In contrast, the QMC calculation irj-tliej 
d ^ oo limit gave a stable AF state for large U > 4tB3'c2l 
We treat here the range of large U {W — 8t) as a 
test case for our analytic method. The calculations were 
performed at a low temperature T — 0.05t which allows 
to describe the magnetic excitations (T ^ J = At'^/U). 
They gave an AF ground state at n = 1 which repro- 
duces correctly the localization of el ectr ons in the limit 
of ?7 —^oo. The magnetization mo (2.2) is only slightly 
reduced by the dynamical effects with respect to its HF 
value, and approaches the HF limit at [/ — > cx3. The 
ground, state is the Neel AF state, as found in the d ^ oo 
limit.c3 Thus, we reduced the self-consiste ntly obtained 
values of the mean-field magnetization mo ( ^.2|) by a fac- 
tor 0.606 in order to simulate the known reduction_nf mo 
by intersite quantum fluctuations in a 2D latticej£3 



m = 0.606mo. 



(4.1) 



After this reduction the calculated values of m approach 
the value of 0.606 in the limit U —i- oo [Fig. ^a)]. |-Qne 
finds also a very good agreement with the QMC dataEl at 
U/t = 2 and 4, and a reasonable agreement at U/t = 8. 




A. Quasiparticles at half filling 

The ground state of the Hubbard model with nearest 
neighbor hopping (i' = = 0) on a square lattice is an 
AF insulator. The insulating behavior and the gap de- 
velop gradually at half-filling with increasing U starting 
from U = Q due to the perfect nesting instability, leading 
to a Slater gap. This gap changes into a Mott-Hubbard 
gap under increasing C/, and the systept. approaches the 
limit of a Heisenberg antiferromagnet.Ell This regime of 
large U was found to be difficult for a-quantitative de- 
scription within the DMFT approchesjij as an accurate 
determination of the energy gains due to AF long-range 



FIG. 1. Antiferromagnetic state for the Hubbard model at 
n = 1: (a) magnetization m ([4.l|), and (b) the AF gap A/J7 
in the 2D Hubbard model, as obtained using the HF approx- 
imation at T = (dashed lines) and DMFT approach at 
T = 0.05t (full lines). The data points in (a) are QMC results 
reproduced from Ref. The diamonds show the results of 
the IPS with the self-energy calculated in SOFT. 

In contrast, the AF gap A is significantly reduced from 
its HF value [Fig. |l|(b)]. This reduction follows from 
a drastic change of the one-particle spectra by dynam- 
ical effects which lead to QP states at the edge of the 
Mott-Hubbard gap which are accompanied by a large in- 
coherent part at higher energies. Also the reduction of 
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the gap found in comes out correctly as shown in 

Ref. For example, we have found a gap of 4.93t at 
U/t = 8, while the corresponding gap in the HF calcula- 
tion is 7.14i. This gap reduction can also be captured by 
the leading dynamical correlations described within the 
SOFT, but only in the regime oi U < 2.5t. The discrep- 
ancy between the SOFT and the DMFT results increases 
with increasing U , with the gap and the magnetization m 
being too small, and finally the AF order disappears and 
the gap closes at ?7 ~ 7t. This shows a very limited ap- 
plicability of the approaches using the self-energy based 
on the SOFTB which are known to underestimate the 
region of stability of magnetic states and fail at large U 
due to the uncontrolled increase of the correlation energy 
in nonmagnetic states. 




-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 

((0-|J,)/t ({J)-|J,)/t ({J)-|J,)/t 

FIG. 2. One-particle excitation spectra as obtained in the 
AF state at n = 1 and T = O.Q5t for the Hubbard model with 
U^8t (t' =t" = 0). 



The spectral functions found within the DMFT (|3^ 
are dominated by the lower Hubbard band (LHB), i.e., 
PES part at w < /i, and the upper Hubbard band (UHB), 
i.e., inverse photoemission (IFES) part at uj > fi, sepa- 
rated by a large gap (Fig. ||). Both the PES and IFES 
spectrum show two distinct energy regimes: (i) narrow 
QF peaks at low energies, i.e., at the edge of the Mott- 
Hubbard gap, and (ii) incoherent and more extended fea- 
tures at higher energies \uj\ > 5t. The overall shape of 
the density of states NLtj^ agrees very well with the ED 
data for a 4 X 4 clusterB The spectra have a character- 
istic k-dependence with the overall weight moving from 
the PES to IPES part along the T ~ X ~ M and F - M 
directions, where F = (0,0), X — (tt, 0) andJVf = (7r,7r), 
in qualitative agreement with QMC dataH The spec- 
tra obey the particle-hole symmetry of the model, with 
spectra symmetric with respect to w = at the X and 
S = (7r/2,7r/2) points. The spectrum at the M point is 
a mirror image of the one at the F point. 

The QF maxima near the Mott-Hubbard gap resem- 
ble those found in the t-J model in ED or within the 
self-consistent Born approximation,El in spite of using a 
local self-energy in the present scheme. This shows that 
the local many-body problem solved within the DMFT 



suffices to capture the low-energy scale relevant for the 
QF propagation. Moreover, unlike in the t — Jz Txm^Ml 
which results in the ladder spectrum for a single holeJ^'E3 
the QF can propagate, as they couple to the spin flips 
of the mean-field bath around site i = at which the 
many-body problem is being solved. The QF dispersion 
is ~ 2 J [Fig. ||(a)], with the maxima along the AF Bril- 
louin zone (BZ^, and remains very close to that found in 
the t-J model.B 



-4.5 




r M X r 

FIG. 3. Momentum dependence in the main directions of 
the 2D BZ, as obtained for the PES spectrum of the Hubbard 
model at half-filling with t' = t" = 0,U = 8t, and T = 0.05t: 
(a) QF dispersion [-E(k) — iJ,]/J; (b) total electron occupation 
(nk) (dashed line) and QF weight a(k) (full line). 

In the HF approximation, the electron occupation fac- 
tors (rik) are larger for the states which belong to the AF 
BZ than for the remaining states outside the AF zone. 
On comparing the weights of the electronic states with 
momenta k and k -I- Q, one finds that also in the DMFT 
the electron weights are much larger within than outside 
of the AF BZ [Fig. |(b)]. The overall PES weight is 
smoothly distributed in the 2D BZ, with the maximum 
(minimum) at the F (M) point, respectively. This result 
agrees well with a QMC simulations, and the present 
data show the same step-like behavior of the electron oc- 
cupation factor (rik) when crossing the X point along the 
F T-T^- M direction as the QMC data at C/ = 4i and 

TABLE I. Values of the model parameters used for 
the presented calculations; the parameter sets chosen 
for La2-i:Sra;Cu04 and YBagCuaOe+i: follow from the 
down-folding procedure of Ref. |^. 



model parameters 


t'/t 


t"/t 


U/t 


J/t 


Hubbard model 


0.0 


0.0 


8.0 


0.50 


La2_a:Sra,Cu04 


-0.11 


0.04 


10.0 


0.40 


YBa2Cu306+. 


-0.28 


0.18 


12.0 


0.33 
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A similar step-like behavior is found as well in the QP 
weight ak along the sa me lin e, determined by integrating 
the spectral functions (3.9) in an energy window of 2J 
which exhausts the range of the QP band in the density 
of states Niuo). The k-dependencc of the QP weight is 
more complex than that of (rik) as two competing effects 
contribute along the T — M and V — X direction when the 
Mott-Hubbard gap is approached: (i) the QP pole moves 
to lower energies and thus the weight increases; (ii) the 
overall PES weight is largest at the F point and gradu- 
ally decreases coming closer to the AF BZ. Therefore, the 
maxima in the QP weight are found close to k = (7r/2, 0) 
and between the F and S = [n 12,11/2) point, while the 
(identical) weights at the X and S point are lower. The 
lowest QP weight is found at the M point, but here in- 
stead a distinct QP exists in the IPES part, in agreement 
with the ED results.Ell 




FIG. 4. Quasiparticles in the AF state at n = 1: the mini- 
mum of the QP band Emin/t (upper part), and the QP weight 
a(k) at k = {n/2, n/2) and averaged over the BZ (lower part) 
as functions of J/t. Filled and empty symbols stand for a(k) 
found in the present DMFT approach and in SCBA of Ref. |5q . 
The inset in the upper part shows Emin/t for < J/i < 10; 
the value of J/t at which the AF order vanishes is indicated 
by an arrow. 

The QP weights a(k) increase with increasing J/t and 
agree surprisingly well with the self-consistent Born ap- 
proximation and ED data for the t-J model in the regime 
of J/t < 0.7, as shown in Fig. ^. The average weight first 
increases somewhat faster than the numerical results of 
Ref. |5^, but then flattens out above J/t ~ 0.6, and sat- 
urates indicating that the t-J model does not represent 
faithfully the hole dynamics in the Hubbard model at 
larger values of J/t, where the excitations to the UHB 
become important. An equally good agreement between 



the self-consistent Born approximation and ED data and 
the present DMFT approach is found at individual k- 
points; the values of o(7r/2,7r/2) are shown in Fig. 
while a very good agreement with ED data at X point 
was presented earlier in Ref. 

The energy at the minimum of the polaron band found 
at the S'-point follows_±he power-law behavior found by 
Martinez and HorschEj in the range of J/t < 0.4 (Fig. 



£^min(kc 

t 



-3.20 



2.94 li 



0.702 



(4.2) 



This power law supports the string picture, but is again 
closer to the full single-hole problem in the t-J model, 
where the data obtained from finite cluster diagonaliza- 
tion could be fitted to the relation -Emin/i = —3.17-1- 
2.93(J/t)°-^3 (ggg j^gf_ tj^an to the t - J, model, 

which gives instead £'min/i = -2^3 -f 2.74( J^/i)^/^ (see 
Ref. ^). It is also quite close to the exact solution of 
the t-J model in the infinite-dimensional lattice, given 
by Ejnin/t = -4Tf2.34(J/i)^/^, which interpolates to the 
Nagaoka state£3 

Finally, we comment on the modifications of the 
spectra introduced by the changes in the parameters 
U and tij. Realistic parameters for La2-a:Sr3:Cu04 
and YBa2Cu3 06+x were estimated both using tl»e cell 
method in the multiband charge-transfer modeljj and 
the down-fokiing procedure in the electronic structure 
calculations .0 Here we use the latter parameters as given 
in Table |, but the sets do not differ significantly. By 
increasing the value of U , one comes closer to the limit 
of the Heisenberg model, and therefore the momentum 
density (rik) is more uniformly distributed over the BZ 
[Fig. |5|(b)]. This quantity depends mainly on the ratio 
of U/t, and thus a similar result is obtained at the same 
value of U with t' and t" non-zero. 
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FIG. 5. Momentum dependence in the main directions of 
the 2D BZ, as obtained for the PES spectrum of the Hub- 
bard model at half-filling with extended hopping parameters 
t' = -0.28t, t" = O.lSt, U = 12t, and T = 0.05t: (a) QP 
dispersion [-E(k) — /i]/J; (b) total electron occupation (rik) 
(dashed line) and QP weight a(k) (full line). 

In contrast, the earlier studies of the t — t' — t" — J 
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model have shown that the dispersion of QP's at low- 
energy are strongly dependent on ther-palues of the ex- 
tended hopping parameters, t' and t" B This strong de- 
pendence is also found in the present calculations based 
on the DMFT approach; the QP's at the S and X points 
are not degenerate any more as soon as t' ^ 0. Here 
we present only the representative result for larger val- 
ues of t' = -0.28t and t" = O.lSt found in SraCuOaCla, 
with minima located close to the X point [Fig. ||(a)]. 
Although the QP weight is dominated by the same com- 
petition between the overall PES weight (rik) and the 
position of the QP maximum with respect to the Fermi 
level, the consequences of sizable t' = — 0.28t are clearly 
visible: the QP weight at the X-point is reduced, and 
the degeneracy of the QP energies found before along the 
T — M and X — T direction [Fig. ||(b)], respectively, is 
now removed. As before, the lowest QP weight is found 
at the M point, and a distinct QP exists in the IPES 
part. Unlike at t' — t" — 0, the latter IPES spectrum 
is different from the PES spectrum at the P point since 
there is no particle-hole symmetry at finite t'. 



B. Spectral properties in spin-spiral states 



As suggested by earlier studies, IllThZI hole doping away 
from half-filling leads to incommensurate magnetic or- 
der. We found the same sequence of spiral phases 
with increasiMg-, doping as in the HF and slave-boson 
calculations the AF order changes first into the SS 
with Q ~ [7r(l ± 2ri),7r{l ± 277)] along the (1,1) direction 
[SS(1,1) state], and then at higher doping into the SS 
with Q = [7r(l ± 277), tt] along the (1,0) direction [SS(1,0) 
state, or an equivalent SS(0,1) state]. The SS states with 
the components of the characteristic Q-vector shifted by 
±277 are physically equivalent and have the same energy. 




FIG. 6. Chemical potential ^/t as a function of doping 5, 
as obtained at T = 0.05t for three sets of parameters given 
in Table |l[ the Hubbard model with U = St (full line) and 
for the model parameters of La2-i:Sra:Cu04 (dotted line), and 
YBa2Cu306+i, (long-dashed line). The regions of phase sepa- 
ration obtained from the Maxwell construction are indicated 
by dashed lines. The inset shows the free energy F/t per site. 



At fixed doping S one finds, however, these phase tran- 
sitions at larger values of U in the present approach which 
includes local correlation-effects, than in the effective 
single-particle theoriesO'Ej This change of the phase di- 
agram follows from t he co rrelation effects which screen 
the value of U to U ( 2.25 ), and strongly depend on the 
magnetic order (Table ||). The highest value of effective 
U/U is obtained in the AF state at half-filling, where the 
double occupancy is strongly reduced and the screening 
is thus ineffective. The screening is stronger in the doped 
cases indicating that the moving electrons correlate and 
avoid each other, leading to much weaker effective re- 
pulsion, and is particularly pronounced in paramagnetic 
states. We found here a surprisingly good agreement for 
the effective interaction U = 1.98t found a,t U — AtMith 
the fitted value of C7 = 2t in the QMC calculations.Ea 

Two regions of phase j^paration which follow from 
the Maxwell constructionliJ were found for the Hubbard 
model at U/t = 8 {t' — t" = 0): a crossover regime 
from the AF to SS(1,1) state for < (5 < 0.11, and 
from the SS(1,1) to SS(1,0) state for 0.22 < S < 0.25, 
respectively (Fig. ||). The value of the chemical po- 
tential fi is U/2 at half-filling, and drops abruptly at in- 
finitesimal doping when it enters the LHB in a doped sys- 
tem. The doping dependence of the free energy indicates 
a phase separation at low doping; this region becomes 
gradually narrower W'th increasing U, in agreement with 
other calculations .EJO In contrast, the transition to the 
SS(1,0) state moves to larger doping with increasing U, 
and finally disappears. Already at the model parameters 
of doped La2Cu04 we found no region of stable SS(1,0) 
state. It is worth noting, however, that in this case a 
small region of almost fiat chemical potential /i was found 
for 5 ~ 1 /8 which could be considered as a precursor ef- 
fect for the phase separation. It might lead P^r^ differ- 
ent magnetic state at still lower temperaturesJl3E3 as the 
stripe structures observed in the neutron experiments. □ 

TABLE II. Values of m agnet ization mo (^^) and the 
renormalized interaction U (2.25), as obtained for the Hub- 
bard model {t' = t" = 0) at T = O.OSi, 5 = I-ti, for different 
magnetic states: antiferromagnetic (AF), spin spiral [SS(1,1) 
and SS(1,0)], and paramagnetic (PM) state. 



ground state 


5 


U/t 


mo 


U/U 


AF 


0.0 


8 


0.871 


0.899 


AF 


0.125 


8 


0.689 


0.755 


SS(1,1) 


0.125 


8 


0.675 


0.735 


SS(1,0) 


0.125 


8 


0.657 


0.733 


AF 


0.250 


8 


0.390 


0.491 


SS(1,1) 


0.250 


8 


0.571 


0.614 


SS(1,0) 


0.250 


8 


0.525 


0.589 


PM 


0.125 


8 


0.0 


0.327 


PM 


0.125 


4 


0.0 


0.494 
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FIG. 7. Spectral functions j4(k,ct;) in the main BZ direc- 
tions: r~X, M~X, and r-M in SS(1,1) state at 5 = 0.125 
and U = 8t for: T = 0.05t (top) and T = 0.5t (bottom). Tiie 
spectra along the T — M direction have been averaged over 
the (1,1) and (1,1) spirals, defined by Q = 7r(l — 2j7)(1, 1) and 
Q = 7r(l + 2r;)(l, 1), respectively. A shadow band below /i in 
A/ — X direction at T = 0.05t is indicated by arrows. 



The k-resolved spectral functions (Figs. ^and|^) allow 
us to identify the generic features of the doped antifcrro- 
magnets described by the Hubbard model, in the regime 
of large U. First of all, the spectra are still dominated 
by the large Mott-Hubbard gap which separates the LHB 
from the UHB. The Mott-Hubbard gap develops from the 
respective gap at half-filling and is considerably reduced 
from U by the QP subbands which occur next to the 
large gap both in the LHB and in the UHB. This large 
gap is accompanied by a smaller pseudogap ~ 2t between 
the occupied (PES) and unoccupied (IPES) part of the 
LHB at low temperature T = 0.05t (taking t ~ 0.4 eV 
it corresponds to ~ 200 K) . This pseudogap results from 
the SS order, and separates the majority and minority 
spin states (with respect to the local coordinates at each 
site) . It is best visible along the T — X and X — M direc- 
tion at 5 — 0.125, and becomes somewhat wider and less 
distinct in the SS(1,0) spiral at higher doping S = 0.25. 
We emphasize that the two features below and above the 
chemical potential originate from the same QP peak at 
half-filling. This shows that the QP found in the spec- 
tral function of one hole in the t-J (or Hubbard) model 
cannot describe the regime with finite doping as the rigid 
band picture breaks down. 




-5 5 10 -5 5 10 -5 5 10 

(co-|i)/t (co-n)/t (co-n)/t 

FIG. 8. Spectral functions A(k, co) as in Fig. but for 
S — 0.25 and (1,0) spiral at T = 0.05t. The conventions are 
the same as in Fig. 

The pseudogap is visible along the T — X direction 
starting from k = (7r/2,0), and the maximum above 
the chemical potential ^ grows gradually towards the X 
point. Consider first the case of lower doping S = 0.125 
(underdoped case). One finds that most of the spec- 
tral weight at the X point is still at w < ^, with a 
sharp QP peak at w ~ —OAAt. Increasing k along the 
X — M direction gives a transfer of the overall weight 
to higher energies, and the QP peak below /i gradually 
looses the intensity, while the peak above fi takes over 
aroimd k — (tt, 7r/3). However, the feature at lo < fi is 
still well visible as a 'shadow' QP band (Fig. 0), with a 
width ~ 2 J. Thus, the QP band of the t-J model is dras- 
tically modified at finite doping, and a new energy scale 
~ 1.5i due to the pseudogap accompanies the dispersive 
QP feature below the chemical potential. 

A similar situation is found also at higher doping 

5 = 0.25 (overdoped regime), and the pseudogap is quite 
pronounced along the T — X and X — M directions (Fig. 
g). However, except for the neighborhood of the F— 
point, more spectral weight is found at high energies. 
Already at the X-point one finds that the peak at cj ~ t 
has a higher intensity than the one below /i. It becomes 
gradually weaker when the M point is approched, and 
disperses in the energy range ~ 3J, while the feature 
below /i still has a similar dispersion ^ 2J as in the 

6 = 0.125 case. We note that the pseudogap increases 
to ^ 2.5t. Moreover, one finds that the QP dispersion is 
broader at 6 = 0.25, indicating the gradual weakening of 
the local magnetic order with increasing doping. 
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The spectra are drastically changed, in particular in 
the low-energy range of \uj — fj,\ < 2t when the tempera- 
ture is increased. At T ~ 0.3t the SS order is unstable 
against the AF order which we interpret as a crossover 
to the small regions of the short-range order with the 
preferably AF ordering of nearest-neighbor spins. The 
spectra found for doping 6 — 0.125 at T = 0.5t con- 
sist of broad maxima which correspond to the LHB and 
UHB, respectively, and only a single maximum is found 
in A(k, uj) next to the X-point. These data, and .also the 
spectral functions for T = 0.33t reported earlier ,c3 agree 
remarkably well with the results of QMC calculations .e3 
The spectra at S = 0.25 and T — O.bt are quite similar 
to those at lower doping S = 0.125, with more weight in 
the IPES part of the LHB. 



3 
< 



N 



X 





-4 -2 
((0-^)/t 



FIG. 9. Spectral functions A(k, u) along two main direc- 
tions in a 2D BZ: F-X [with the step of (7r/3,0)], &ndX-Z, 
where Z = (7r,7r/2) [with the step of (0,7r/6)], in SS(1,1) state 
as obtained for the model parameters of doped La2-i:Sri:Cu04 
(Table |) at doping S = 0.125 (left) and S = 0.25 (right), and 
after averaging over all equivalent SS states with different val- 
ues of Q. The dispersive feature with the strongest intensity 
is indicated by vertical lines. 

We do not intend to present detailed analysis of the 
spectra obtained using the extended hopping param- 
eters which correspond to the electronic structure of 
La2_a;Sra;Cu04 and YBa2Cu306+a:, respectively. In- 
stead, we point out the important similarities and dif- 
ferences to the Hubbard model as far as the SS states 
are concerned. Consider first the effective parameters of 
La2-a;Sra;Cu04. First of all, a narrow QP band is also 
found below the Fermi energy (Fig. |), but the mea- 
sured dispersion between the F and X point is ^ 0.80t 
(- 2.56t) at 5 = 0.125 {S = 0.25), while it amounts only 
to ~ 0.36t at S — 0. Note that the energies of the QP 
peak are much closer to each other for the present pa- 
rameters than in the Hubbard model with t' = t" = 0, 



where one finds instead the dispersion of 1.37t (3.58i) 
at 5 = 0.125 {5 = 0.25) at U = 8t, while it amounts 
to Lit in the undoped case. This gradual widening of 
the QP dispersion with increasing doping may be under- 
stood as a consequence of the admixture of ferromagnetic 
components with increasing doping in the SS(1,1) states. 
The same trend is also observed for the parameters of 
YBa2Cu306-i-a;, where one finds the QP states in the PES 
part separated by ~ 0.39t 1.94i) between the F and 
X point at 5 = 0.10 (0.25) (Fig. while this splitting 
is only ~ 0.07< at half-filling. 

Finally, the finite hopping elements to the second and 
third neighbors stabilize the SS(1,1) state with respect 
to the SS(1,0) state also at higher doping (see Fig. ||), 
and therefore the intensity at the X point does not cross 
the Fermi level even at 6 = 0.25 for both parameter sets. 
In fact, taking J = 0.125 meV {t/J = 3), the QP state 
at the X point is found at a; ~ —0.56 eV, and does not 
change significantly as a function of doping (Fig. |o|). 
In contrast, in the ARPES experiments for Bi2Sr2CaCu2 
the QP state at X-point is found at energy ~ —0.20 
eV (~ — O.QSjB eV) in the underdoped (optimally doped) 
compound.E3 This indicates that either an improved so- 
lution of the many-body problem is still required, or the 
actual magnetic order in these compounds might be dif- 
ferent from SS states. However, the observed increase of 
the onset of incommensurability with increasing U and 
t' is consistent with the observations made by Igarashi 
and FuldeEII and with QMC calculations of Duffy and 
MoreoS 
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FIG. 10. The same as in Fig. but for the model pa- 
rameters of doped Y(Bi)-superconductors (Table |) at doping 
S — 0.10 (left) and 5 — 0.25 (right), and after averaging over 
all equivalent SS states with different values of Q. 
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C. Total densities of states 



We already pointed ou a very go od ag reement be- 
tween the calculated densityi-Qf states ( 3.11 ) and the re- 
sults of ED by Dagotto et a/.E3 Here we present instead a 
comparison between the density of states obtained within 
the DMFT method and that found in the HF approxi- 
mation (Fig. 



11) 



First of all, one notices a narrower 
gap which separates the QP subbands in the DMFT of 
width ^ 2 J, instead of the HF one-particle states, on the 
scale of ~ 2t. This part of the spectral density might also 
be reproduced in effective single-particle approaches, as 
for instance in the slave-boson mean-field theory. How- 
ever, the incoherent parts which extend on the energy 
scale down to |w — /u| ~ 9t result from many-body scat- 
tering and can only be reproduced if the dynamical part 
of the self-energy is included. The overall width of the 
subbands at w < and w > is ~ 7t, respectively, as 
known from the analysis of the t-J model in ED and in 
QMC calculations Jj 

It is evident that due to the changes of N{uj) in the 
range of — /i| < 1.5t with respect to QP band in the 
undoped system, the low-energy part of the spectrum 
cannot be reproduced in a renormalized one-particle the- 
ory. The pseudogap in the doped systems is not visible 
in the HF density of states, and it remains a challenge 
whether an effective one-particle theory which captures 
this essential new energy scale could be constructed. As 
expected, the agreement between the HF and DMFT 
density of states improves somewhat at higher doping 
S = 0.25, where the Mott-Hubbard gap is gradually lost, 
and the system approaches the single-particle limit. We 
note, however, that the gap between the LHB and the 
UHB relies in our approach on the magnetic order, and 
a more accurate approach in the strongly doped regime 
at large U would instead have to include the scattering 
on local moments. 

In spite of a very good agreement_between the present 
DMFT approach and the ED dataj^j it is interesting to 
investigate to wh at extent the analytic formula for the 
self-energy ( 2.22 ) describes the hole dynamics in a doped 
system. Therefore, we performed also a DMFT-QMC 
calculation of local 'E,Q{iuj) for SS states, and the corre- 
sponding densities of states, shown in Fig. ^ The QP 
peaks are very close to each other at half-filling, while 
the incoherent states at higher energies in the LHB and 
UHB have almost the same weights, but are moved to 
somewhat higher energies in the QMC calculation. The 
increase of the spectral weight close to the Fermi level is 
well pronounced in the latter calculation at 5 = 0.125, 
but one finds instead a pseudogap smaller by a factor 
close to five. 
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FIG. 11. Total densities of states N{uj) as obtained within 
DMFT (full lines) for 5 = (AF state), 0.125 [(1,1) spiral], 
and 0.25 [(1,0) spiral] with U/t = 8 and T = 0.05t. The 
dashed lines show N{uj) for the magnetic ground states found 
in the HF approximation. 



However, one should realize that the present calcula- 
tion performed at low temperature T = 0.05i corresponds 
in practice to the ground state, while the same tempera- 
ture in QMC includes already thermal fluctuations which 
considerably reduce the size of the pseudogap. Indeed, 
we find using the ED method to solve the self-energy 
within DMFT the pseudogap in the SS state of ~ 0.7t. It 
might be expected that this reduction of the energy scale 
would result in a better quantitative description of the 
spectral functions and the related excitations across the 
pseudogap, leading to a reduced energy scale for the low- 
energy features of the optical conductivity (Sec. 0). We 
also found more extended energy range of the incoherent 
states which belong to the UHB in the QMC calculation. 
Altogether, the comparison with the DMFT-QMC calcu- 
lation demonstrates that the analytic method developed 
in this paper is very useful for a qualitative insight into 
the possible changes of magnetic states under doping and 
their consequences for the properties measured in exper- 
iment. 
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FIG. 12. Total densities of states N{ij) as obtained within 
DMFT at ?7 = 8i by calcu lating the self-energy either using 
an analytic formula (^.22|) (solid lines), or by DMFT-QMC 



method (dashed lines). Different panels show the results ob- 
tained for 5 = (AF state), 0.125 [(1,1) spiral], and 0.25 [(1,0) 
spiral], respectively, at T = 0.05t. 



V. OPTICAL AND TRANSPORT PROPERTIES 

A. Optical conductivity in the Hubbard model 

The evolution of the spectral functions ^(k, uj) and 
the d ensit y of s tates N{uj) with doping, reported in Sees. 
[V A and III A , motivates an investigation of the optical 



properties. Here we make use of the theory introduced 
in Sec. IIIB, where we have shown how the optical con- 
ductivity can be derived from the local self-energy in the 
present DMFT treatment. 

As an illustrative example, we concentrate on the op- 
tical conductivity found for the Hubbard model with the 
nearest-neighbor hopping {t' — t" — 0) and U ~ 8t. We 
present the optical data in Figs. ^ and |lj at two tem- 
peratures: T = O.OSi and 0.2t. While the magnetic order 
is AF at half-filling, the SS states characterized by the 
wave-vector Q = [7r(l ± 2rix),n{l ± 2riy)] change with 
doping and temperature. At lower doping S = 0.125 we 
find a SS(1,1) with r]^ = r/y ^ 0.125 (0.09) at T = 0.05i 
(T = 0.2t), respectively, while at higher doping S = 0.25, 
a SS(1,0) state {r]y = 0) with r]^ = 0.25 (0.23), or an 
equivalent SS(0,1) state, is found instead. 
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FIG. 13. Optical properties as functions of energy uj/t for 
the Hubbard model with i7 = 8t at low temperature T = 0.05t 
for (5 = (dashed lines), 5 = 0.125 ( full lines), and 5 = 0.25 
(dashed-dotted lines): (a) real part of the optical conductivity 
<j'{i-o); (b) imaginary part of the optical conductivity a"{Lo); 
(c) scattering rate l/r(aj); (d) effective mass m* [ui) /nie. 



At half- filling one finds a large gap below ~ 4.9i at 
U /t = 8 and no Drude peak which shows that the system 
is in the insulating phase.E^ The conductivity at w > 4.9i 
is incoherent and originates from the excitations across 
the Mott-Hubbard gap. This changes drastically when 
the system is doped and two new features occur at lower 
energy: the Drude peak, and the mid-gap state at w ~ 2t, 
both with the increasing intensity between i5 = 0.125 and 
5 — 0.25 at low temperature (Fig. |l^). These features are 
accompanied by an incoherent background of the excita- 
tions within the LHB. The peak at ~ 2t corresponds 
to excitations across the pseudogap; as such it is more 
influenced by the increasing temperature in the under- 
doped regime, where the SS(1,1) state is less robust than 
the SS(1,0) state in the overdoped regime. 



Below Lo = AM the frequency dependent scattering 
rate 1/t(ijj) and the effective mass m*{uj)/me can also be 
divided into two regions: (i) above ~ 2i the scattering 
rate increases monotonically with increasing frequency, 
(ii) below w ~ 2i it has a maximum at energy lu ~ l.l5t 
(IM) for S = 0.125 (0.25), and drops to zero for w ^ 
at finite doping. This behavior for uj ^ and T ^ is 
consistent with the Fermi liquid behavior whic h foll ows 
from the local approximation to the self-energy (2.11). A 
finite value at a; = is a numerical effect due to finite 
broadening of the spectra (e = O.lt). 
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FIG. 14. Optical properties as functions of energy ijjjt for 
the Hubbard model with U = 8t at intermediate temperature 
T = 0.2t; the meaning of lines as in Fig. 

The frequency region in which the scattering is sup- 
pressed has a direct relation to the existence of a pseudo- 
gap region in the singl e-particle spectral function A(k, iS), 
reported in Sec. IV B , and indicates that SS LRO reduces 
the scattering of the charged carriers in the energy range 
uj < 1.15t. At the same time, the effective mass m*{u}) 
rises to a maximum value of ~ brue within the pseudogap 
region, and is found to be rather independent of hole dop- 
ing. As the temperature increases to T = 0.2t, the pseu- 
dogap disappears and the region of suppressed scattering 
is filled up in the underdoped regime with S = 0.125, 
while the low scattering persists for lu < l.Ot at 6 = 0.25 
(see Fig. At the same time, the mid-gap state in 

the real part of the optical conductivity changes into a 
smooth feature which extends down to the Drude peak 
for 5 = 0.125, contrary to the case with 6 = 0.25 where 
the spectral weights of the above two features remain 
well separated. This is clearly related to the behavior 
observed in A(k, lo) with increasing temperature, where 
the pseudogap along the X — M direction filled up with 
spectral weight as T increased for 5 = 0.125. 

These changes of the mid-gap state with temperature 
are due to the changes of the magnetic correlations in the 
doped systems included in our calculations which do not 
distinguish between long-range and short-range magnetic 
order, but treat local dynamical correlations. However, 
there are indications that the mid-gap feature results 
from an interplay betweea short-range magnetic order 
and electron correlations .Eil Therefore, the rather strong 
evolution of the low-energy weight with increasing tem- 
perature shown in Figs. |l^ and |l^ may be overestimated 
in the present treatment of the self-energy which does 
not allow to get a metal-insulator transition without an 
accompanying magnetic LRO. We also note that the mid- 



gap states are likely a bare consequence of the strongly 
correlated nature of oj^tical and one-particle excitations 
in the Hubbard modelj^ and it is still a challenge to de- 
scribe them better in a theory which would treat the AF 
and paramagnetic states with local moments on equal 
footing. 

The frequency-dependent scattering rate allows us to 
find a crossover temperature T* at which the pseudogap 
closes. We estimated that T* ~ 0.26i for S = 0.125, 
and observed a monotonic increase of 1/t{uj,T*) up to 
Lo ^ A.lt. At T = 0.2i the effective mass increases up 
to ~ IOtoe within the pseudogap at (5 = 0.125. At half- 
filling and T = 0.2i one finds that the charge-transfer 
gap is only slightly reduced from its value at T = 0.05^, 
and the insulating behavior is accompanied by AF LRO. 
We estimated the Neel temperature for U = 8t to be 
Tjv ~ 0.62 J. 

Further evidence for a characteristic crossover temper- 
ature T* m ay be found in the behavior of the in-plane dc 
resistivity (3.27). The resistivity received a lot of atten- 
tion in connection with the observed normaL-state pseu- 
dogap in the electronic excitajtion spectrumJS and from 
theoretical point of view£3'tJ In fact, the physical ori- 
gin of linear T-dependcncc of p{T) for samples of high- 
Tc compounds close to the optimal doping level remains 
puzzling. 




FIG. 15. Resistivity p(T)/po as a function of temperature 
T/t as obtained for the Hubbard model with U = 8t for 
S = 0.125 (full line) and S = 0.25 (dashed line). The in- 
set shows the weight z*{T) ( |5.l[ ) found at the X point at 
S = 0.125 (full line), and averaged weights over the Bril- 
louin zone at 5 = 0.125 (dashed line) and at 5 — 0.25 
(dashed-dotted line). Arrows in the inset indicate T* . 

The results for p^x (T) obtained for the Hubbard model 
at two doping levels, 6 = 0.125 and S = 0.25, are shown in 
Fig. |l^. At low temperatures, T < 0.06t, the resistivity 
shows Fermi-liquid behavior for both hole densities, i.e_ 
Pxx{T) (X T^. As usually in the DMFT calculations ,EI 
the T^-dependence of pxx{T) originates from the low- 
frequency behavior of the imaginary part of the local self- 
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energy. In the regime of high temperatures (T > 0.9<), 
the resistivity increases hnearly with temperature which 
is due to temperature independence of the spectral func- 
tions Ao-o-' (k, cj) at T — *■ oo, and the high temperature 
hmit of the derivative, {~dnp{u!)/duj) 1/(4T), thus 
leading to (Txx{T) oc 1/T, i.e., Pxx{T) oc T. As the 
temperature is lowered, the magnetic moments gradually 
build up, and a kink in the resistivity appears. There- 
fore, the increase in the resistivity as the temperature 
is lowered can be attributed to the enhancement in the 
scattering of electrons by local spin fluctuations. 

On the contrary, for large hole doping the system is 
a better metal, hole spin correlations are gradually lost, 
and the increase of resistivity is less pronounced in the 
temperature region T < 0.88i at 6 = 0.25. The maxi- 
mum of Pxx{T) for 5 = 0.125 is located almost exactly 
at T w 0.26t 750if taking the experimental value of 
the superexchange J = 125 meV), where the pseudogap 
in the single-particle excitation spectra opens leading to 
a suppression of the effective scattering rate 1/t(ijj,T), 
as discussed previously. This defines the crossover tem- 
perature T* . Remarkably, the change from a linear to a 
nearly-hnear T-dependence, Pab{T) oc T^+' (e > 0), of 
the in-plane dc resistivity of La2-2;Sr2;Cu04 was found 
to be at T* ~ 600 K for a; ~ 0.13 and was attributed 
to the openina-pf a pseudogap in the electronic excita- 
tion spectrumEa However, the saturation of the resistiv- 
ity Pab{T) cannot be observed in a real system as the 
carriers also couple to other bosonic excitations, e.g. to 
phonons, which are neglected here. 

Upon lowering the temperature below T*, at T < 0.24i 
for both 5 = 0.125 and for 5 = 0.25, we observed a nearly 
linear T-dependence of Pxx{T). In this temperature 
range the SS wave vector, Q = [7r(l ± 2772,), 7r(l ± 2?7y)], 
becomes strongly temperature dependent, and main- 
tains the directional deviation from the AF wave vector, 
Qaf = (tt, tt), with rjx = riy = rjiT) for 5 = 0.125 and 
ry^ = r]{T) {riy = 0) for 5 = 0.25. In both cases t]{T) in- 
creases from ri{T*) ~ with decreasing temperature and 
saturates at its ground-state value rjiT = 0) ~ J below 
T ~ 0.08<. In the hnear regime (T < T*) the resistivity 
can be fitted quite well laz^a linear T-dependence, as ex- 
pected for the SS statesjlil p^* (T) = p^* (0) + (poS'^T, 
with p^i{0)/po = -1.05 (-0.25) for S = 0.125 (0.25), 
respectively, where the increase of the negative temper- 
ature coefficient /c^^(O) is a further manifestation of the 
gradual loss of local magnetic moments as doping is in- 
creased. On the contrary, in the paramagnetic phase 
o£_the Hubbard model at d = cxd one finds p^^(O) > 
O.lZI Furthermore, the slope of Pxx{T) in the low tem- 
perature regime is given by ^ ~ 1.46 independent of 
hole density. This value is larger by about a factor of 
2.5 than the respectpp slope found in the retraceable- 
path approxiHjation,Lj and in the ED studies at finite 
temperature)^ being = 0.55 and 0.60, respectively, and 
demonstrates that the changes in the magnetic order with 
increasing temperature influence significantly the system 
resistivity. Unfortunately, such effects cannot be studied 



in the ED method due to the small size of considered 
clusters. 

In order to further support our observation that the 
crossover temperature T* is related to the pseudogap in 
the single-particle excitation spectrum we plot in the in- 
set of Fig. ^ an average of the single-particle spectral 
weight within an energy pstipdow cx T around the Fermi 
energy a; = 0, defined by,OL3 



^(T) = -^C-Q(k;f,r = /?/2) 
A(kF,w) 



duj 



cosh(/3w/2)' 



(5.1) 



Similarly, a measure for the temperature dependence 
of the density of states at the Fermi energy Njfi) 
is obtained from the local Green's function ( 2.11 ), 
^IciT) = -^EaG.aQir = (3/2). In the low tem- 
perature limit A^(0)|-C|an be obtained from the relation 
iV(0) ~ (3zl^iT)/njB which gives ~ 0.20 (~ 0.26) for 
6 = 0.125 {S = 0.25), respectively. However, one finds 
that the one-particle density of states at the Fermi en- 
ergy does not evolve smoothly to the low temperature 
values, but instead states are depleted from the region 
u; ~ p as T is reduced below T* ~ 0.26i (0.43t) for 
6 = 0.125 {S = 0.25), respectively. In particular, we ob- 
served a faster loss of the QP weight with momentum 
kx = (7r,0) for S = 0.125 (Fig. This shows that the 
opening of the pseudogap in the one-particle excitation 
spectrum at (vr, 0) coincides with the suppression of the 
effective scattering rate l/r(ti;,T). 

Experimentally, the resistivity changes from a linear to 
a nearly-linear T-dependence at T* of the order of 500 
K. Although our calculations do not allow to interpret 
the linear part of Pxx{T) at high temperature T > T* 
as only the electronic degrees of freedom are included, 
we note that the enhanced slope of Pxx(T) at low tem- 
perature T ^ 100 K and the negative temperature coef- 
ficient agree qualitatively with the experimeatal results 
for YBa2Cu307_a; in the underdoped regime.Ej Our cal- 
culations confirm the conjecture of Shraiman and Siggia 
of a nearly-linear T-dependence_of the resistivity for a 
system with SS magnetic order .lj These features can be 
seen as generic fingerprints of incommensurate magnetic 
correlations. 
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B. Implications of extended hopping 
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FIG. 16. Optical properties as functions of energy u/t for 
the Hubbard model with extended hopping parameters of 
La2-a:Sri:Cu04 (Table |) at low temperature T = 0.05t for 
5 = (dashed lines), 5 = 0.125 ( full lines), and 5 = 0.25 
(dashed-dotted lines): (a) real part of the optical conductivity 
cr'(a;); (b) imaginary part of the optical conductivity a"{u>); 
(c) scattering rate l/r(aj); (d) effective mass m* (u) /me. 

Similar changes in the optical excitation spectra as a 
function of hole doping were also found using the effec- 
tive single-band models with the parameters represen- 
tative for La,2-xSTxCu04 (Fig. |l^ and YBa2Cu306+2: 
(Fig. |l^), respectively. Due to somewhat larger values of 
the effective U, the gap in the optical spectra increases 
to ~ 6.5i and 7 At in these two compounds. One finds 
again that the Drude weight and the mid-gap state ap- 
pear in the conductivity of doped systems. For the pa- 
rameters of La2-2:Sra;Cu04 (YBa2Cu306+2:) the region 
of suppressed scattering extends up to ~ IM (~ IM) 
at 5 — 0.125. This regime of low w gives an enhanced 
effective mass ~ 4me for both sets of model parame- 
ters. At larger doping 5 = 0.25 the coherence of the 
charge carriers is enhanced by t' and t" hopping, and 
one finds a significantly reduced effective scattering be- 
tween charged carriers, extending with roughly no struc- 
ture over a rather broad energy range. Simultaneously, 
the effective mass ~ 1.5me is only little enhanced at low 
energies. 

The overall shape of u'{lo) (Fig. |l6|) shows a quali- 
tatively similar behavior to the optical cmiductivity of 
La2-a;Sr2;Cu04 reported by Uchida et aZ..c3 At low dop- 
ing the mid-gap band centered at uj ~ l.lt (correspond- 
ing to 0.53 eV for J = 125 meV and the present param- 
eters with J — QAt) is clearly distinguishable from the 
Drude contribution. It moves to higher energy uj ~ 2.2t 



(0.7 eV) &t 5 = 0.25. It is quite remarkable that our 
DMFT calculations reproduce qualitatively the struc- 
tures observed in the frequency dependent effective scat- 
tering rate l/T(tj)|-and in the effective mass m*(a;)/me 
of La2_a;Sra;Cu04.c3 In particular, the strong doping de- 
pendence of 1 /t{uj) and m* [lj) / me show the same trends, 
namely a pronounced reduction of scattering and effec- 
tive carrier mass for the heavily doped systems, and fur- 
ther justifies the importance of extended hopping param- 
eters in the cuprates. This behavior originates from an 
increase of QP weight in the single-particle excitation 
spectrum induced by doping. 
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FIG. 17. Optical properties as functions of energy uj/t for 
the Hubbard model with extended hopping parameters of 
YBa2Cu306+x (Table % at low temperature T = 0.05f for 
5 = (dashed lines), 5 = 0.10 ( full lines), and S = 0.25 
(dashed-dotted lines): the meaning of different panels is the 
same as in Fig. |l^. 

Puchkov et a/0 reported extensive studies of the in- 
frared properties of YBa2Cu306+a;, Bi2Sr2CaCu208+a;, 
and other high- Tc compounds. They found that the far- 
infrared effective scattering rate 1/t{ijj) and the effec- 
tive mass m*{uj)/me differ significantly between under- 
doped and optimally doped samples above Tc. The opti- 
mally doped samples show a structureless and lower ef- 
fective scattering rate 1/t{lu) and a nearly constant and 
unrenormalized mass r7i*(a;). On the contrary, in the 
underdoped samples the scattering between the charged 
carriers below « 0.12 eV is strongly suppressed and 
TO*(aj)/TOe is enhanced in the low-energy region. These 
observations are in remarkably good agreement with our 
findings and supports our conclusion that the observed 
doping dependence of l/r(w) and m*{uj)/me originate 
from an increased coherence of the one-particle excitation 
spectra reported in [VB| , and experimentally observed in 
ARPES spectra of Bi2Sr2CaCu208+a; by Kim et al.B 
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The suppression of 1/t{uj) below w ~ 0.12 eV originates 
from the opening of a pseudogap in the one-particle ex- 
citation spectrum. Using J = 125 meV and the value 
of J = t/3 adequate for YBa2Cu306+a;, we find the en- 
ergy threshold below which QP scattering is strongly sup- 
pressed in the weakly doped system at ~ 0.68 eV. Un- 
fortunately, this is about a factor of five larger than the 
experimental value for underdoped YBa2Cu306+2; being 
J. Similar discrepancies in the energy of t he QP state 



with momentum (7r,0) were reported in Sec. IVB 



C. Drude weight and spectral weight transfer 



Finally, we compare the Drude weig ht D ( ^.26D and the 
kinetic energy density associated with x-oriented links 
i—kx) ( 3.16| ) for the three model parameter sets in Fig. 
|l|. The DMFT gives (-fc^) = 0.46t for the Hubbard 
model at half-filling {5 — 0) with U = 8t which is by a 
factor ~ 1.81 smaller than the HF result and is in excel- 
lent agreement with the value of {—kx) = 0A9t obtained 
in QMC calculations .Ea We also found an overall satis- 
factory agreement of {—kx)r^ a function of doping with 
ED data of Dagotto et aLE2 The kinetic energy {—kx) 
increases with doping not only because the actual carrier 
density changes, but also as a consequence of changing 
wave vector in the SS state Q = [7r(l ± 2?7),7r(l ± 277)] 
with a gradually increasing pitch 77 allowing for coher- 
ent electronic transpoijt-|through the system. In agree- 
ment with QMC dataj23 we observe that increased ex- 
tended hopping amplitudes accelerate the SS formation 
and result in a stronger increase of {—kx) with 5 for the 
YBa2Cu306-i-a; than for La2-£i;Sra;Cu04 model parame- 
ters. The a;-directed kinetic energy shows a linear doping 
dependence in the regime of low hole doping, and {—kx) 
changes by a factor of ^ 1.66 (1.22) with respect to half- 
filling in the case of the YBa2Cu306+a: (La2-a;Sr2:Cu04) 
model parameters, one finds an faster increase of total 
spectral weight in the case of stronger hopping to second 
and third neighbors, as realized in YBa2Cu306+x- 
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FIG. 18. Kinetic energy along the a;-direction {—kx)/t and 
the Drude weight D/t as functions of the hole doping S for 
representative values of parameters given in Table ^: Hubbard 
model (circles), La2-a;Sra;Cu04 (squares), and YBa2Cu306+a; 
(diamonds). 

The calculated total optical spec tral weights are cx 
{—kx) following the optical sum rule (3.22), and we made 
a quantitative comparison with the experimental data. 
The doping dependence of the total integrated spectral 
weight below the charge-ittansfer band edge at 1.5 eV 
reported by Cooper et alE3 for La2-3;Sr3;Cu04 is strik- 
ingly similar to the numerical data of Fig. 
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model reproduces a rapid increase of spectral weight up 
to ~ 10% Sr-doping and a rather doping independent 
spectral weight in the range of 0.1 < a; < 0.2. The in- 
crease of oc {—kx) with increasing doping is faster for the 
parameters of YBa2Cu306+a;, with the integrated spec- 
tral weight increased by ^ 1.7 at S — 0.18 with respect to 
its value at 5 = 0. This value compares again very well, 
taking the simplicity of the effective single-band Hubbapd 
model, with a factor of ^ 1.8 found by Orenstein et aZ.Ej 
in the compound with highest Tc- 

At (5 = wc find a vanishing Drude weight for all three 
sets of model parameters, and the system is an insulator. 
This is of course an expected result at half-filling, but in 
the present context it serves as a test of the internal con- 
sistency oftheory, as the kinetic energy term {—kx) 7^ 
in Eq. ( |3.26 ), and has to be compensated by the current- 
current correlation function Axxi<i = 0, 2'KiT) in the limit 
of low temperature. At small hole doping we observed an 
almost perfect linear increase of the Drude weight with 
i5 for all three sets of model parameters which is an in- 
dication of strong eleqti-OH correlations near the Mott in- 
sulator at half-fiUing.EJO. Such a behavior is compati- 
ble with a picture of a dilute hole gas in a background 
with SS LRO which contributes to the optical response. 
However, the crossover to a metal due to increasing dop^ 
ing has been analyzed recently using scaling tlif!ory,Ej 
and ED technique combined with scaling theory,L3 which 
give D (X 5^ for small doping concentration (5 in a 2D t-J 
model. This last result is in sharp contrast to the present 
picture of a dilute hole gas in an AF or SS background, 
and might indicate that other correlations are realized in 
the spin background when the system is doped, namely 
that the dilute hole gas is unstable towards microscopic 
phase separation, such as realized in polaronic solutions 
or stripe phases. 

The present results demonstrate a substantial trans- 
fer of spectral weight to low energyj-jn the doped sys- 
tems. We already pointed out earlicro that the spectral 
weight transferred into the LHB in the one-particle spec- 
tra agrees with the predictions of. perturbation theory 
in the strongly correlated regimeO In the optical spec- 
tra for the Hubbard model at U/t = 8 one finds that 
the weight transferred into the region below the Mott- 
Hubbard gap is increased by a factor ~ 1.3 with respect 
to S = 0.125 when the system is doped to 6 = 0.25. 
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This change is significant as the total weight obtained 
from ( 3.22 ) via {—kx) remained roughly constant [see 
Fig. |l8|(a)], indicating a spectral weight transfer from 
the high- to the low-jejjiergy region in the single-particle 
excitation spectrum.o In particular, the weight trans- 
fer is in favor of the Drude weight [Fig. ^(b)], which 
increased in the same doping range by a factor ^ 2.15 
although the hole density increased only by a factor 2. 
These changes in the coherent optical wei ght a re consis- 



tent with the observation made in Sec. IV B that the 
single-particle excitation spectra become more coherent 
as the hole density increases. 



VI. SUMMARY AND CONCLUSIONS 



We reported a generalization of the DMFT to the mag- 
netically ordered states and showed that this method al- 
lows for a very transparent study of spectral properties of 
the Hubbard model at and close to half-filling. The cru- 
cial step is the derped formula for the self-energy using 
the Berk-SchriefferEa spin-fluctuation exchange interac- 
tion with aa effective potential due to particle-particle 
scattering. a We have demonstrated that this treatment 
of the many-body effects reproduces the leading depen- 
dence on doping and temperature and gives a very favor- 
able comparison with the available numerical data ob- 
tained in the QMC and ED calculations for a 2D Hubbard 
model. Although the k-dependence of the self-energy was 
not included, the spectral functions for a single hole in 
a Mott-Hubbard insulator -agrees well with the known 
structure for the t-J model,El and gives the PES spectrum 
consisting of a coherent QP peak with a dispersion ~ 2 J, 
and an incoherent part of width ~ It at lower energies. 
We have verified that the QP weight,-^rees well with the 
ED data in the. range oi J/t < 0.7j^ and supports the 
string picture.LJ Furthermore, the calculation reveals a 
nontrivial relation between the electron occupation fac- 
tors, (fT-k), and the QP weights, Ok, and shows that the 
maximum of ak is shifted away from the (7r/2, 7r/2) point, 
in agreement with the ED results of Eskes and Ederj£3 

Our study has shown that doping of a Mott-Hubbard 
insulator leads to an incommensurate magnetic order at 
low temperatures, which depends on the actual values of 
the hopping parameters and the Coulomb interaction U . 
This kind of magnetic order induces a pseudogap in the 
one-particle spectra which is one of the generic features 
of the doped Mott-Hubbard insulators. The dependence 
of the pseudogap on the incommensurate magnetic or- 
der explains why it could not be obsecved in ED data on 
small clusters at finite temperaturej£j or in the infinity-, 
dimensional Hubbard model in the paramagnetic state.ESI 
This new energy scale due to a pseudogap of magnetic ori- 
gin demonstrates a combination of physics arising from 
the Slater picture and the Mott-Hubbard description of 
strongly correlated electron systems. 

The coherent QP states survive in the doped systems. 



in agreement with the QMC and ED results. However, 
the numerical studies suggest that a strong k-dependence 
of self-energy might be necessary to describe the spectra, 
as the QP dispersions change. This failure of the rigid QP 
band picture finds here quite a different explanation: the 
changes of the QP dispersion follow from the incommen- 
surate magnetic order which develops with doping, and 
the leading effects in the hole dynamics are still captured 
by a local self-energy. 

The one-particle and optical spectra are interrelated, 
and the opening of a pseudogap at low temperatures leads 
to a mid-gap state next to the Drude peak in the optical 
conductivity, both with growing intensity under increas- 
ing doping. Such features, observed in the SS states at 
low temperatures, as the suppressed scattering rate and 
large effective mass in the underdoped regime, and al- 
most no enhancement of the effective mass in a broad en- 
ergy range in overdoped systems are in remarkably good 
qualitative agreement with the experimental findings in 
the cuprates.oO This is consistent with the reduced den- 
sity of states N{fj,) at the Fermi energy at low temper- 
ature. With increasing temperature the value of N{ii) 
increases, which could not be explained in paramagnetic 
calculations performed within the DMFT approach. It 
should be realized that such a strong temperature de- 
pendence of N{p) should have important consequences 
for several measurable quantities in the normal phase, as 
for example Knight shift. 

Here we limited ourselves to the qualitative conse- 
quences of the extended hopping t' and t" for the one- 
particle and optical spectra. First of all, the QP disper- 
sion is strongly influenced by these parameters, and at 
half filling reproduces the experimental psadth and dis- 
persion of the QP band in Sr2Cu02Cl2O Second, the 
deviation of the characteristic Q-vector from the AF vec- 
tor (tt, tt) increases as a function of doping in SS states, 
and this process is accelerated by a finite value of second- 
(i')and third-neighbor (t") hopping. This explains why 
the systems with extended hopping are more metallic 
which is indicated by the low effective mass and larger 
Drude weight. 

The dependence of the magnetic order on tempera- 
ture has also rather drastic consequences for the mea- 
surable quantities. The onset of magnetic order below a 
characteristic temperature results in quite different one- 
particle and optical spectra at low temperatures from 
those obtained in a paramagnetic phase. The changes of 
the spiral Q-vector with decreasing temperature allow to 
introduce a crossover temperature T* , below which the 
low-lying excitations are gradually modified along with 
the changes in local magnetic order. Such a modification 
gives a quasi-linear resistivity, and verifies the conjecture 
of Shraiman and Siggia.U 

In spite of very good agreement for the undoped sys- 
tems, however, we identified several important features 
which do not agree with the experiments in the doped 
cuprates even on a qualitative level, that might indi- 
cate that either a more accurate treatment of the many- 
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body problem is necessary, or more complex magnetic 
structures are stabilized in these compounds: (i) The 
value of the pseudogap in the one-particle spectra and 
the accompanying energy scale for the suppressed scat- 
tering rate in the optical conductivity arc overestimated 
by a factor close to five with respect to the experimen- 
tal observations, (ii) The SS(1,1) state obtained for 
the La2-£!;Sra;Cu04 model parameters leads to a differ- 
ent splitting of the magnetic scattering peak in neutron 
experiments than the experimentally observed (1,0) and 
(0,1) splittings.Ll (iii) The incommensurate order deviates 
too fast from the AF state for the model parameters of 
YBa2Cu306+£c which results in different spin-spin cor- 
relations than those observed in experiment, ancLa QP 
peak at the X-point moving to too low energies.^] (iv) 
The doping behavior of the pseudogap and of the related 
crossover temperature T* is opposite to the one observed 
in the cuprates. In these materials the pseudogap and 
T* decrease upon doping, whereas here the correspond- 
ing quantities increase from the S — 0.125 to the S — 0.25 
case. With increasing doping charge fluctuations become 
more and more important which gives rise to the sup- 
pression of magnetic order, and consequently the pseu- 
dogap closes. However, such correlations are underesti- 
mated in the present treatment, and one finds instead 
a persisting pseudogap. (v) Finally, the spiral spin or- 
dering in the (1,1) direction contrasts with experimental 
evidences from neutron scattering in the cuprates, sug- 
gesting that stripe ordering might play a pmpunent role 
in these systems at very low temperaturesJaH We have 
found a phase separation at low doping levels and there- 
fore the presently studied dilute hole gas in SS states 
is unstable towards magnetic polarons or stripe phases 
at doping levels lower than S ~ 0.1. This motivates a 
further search for more complex magnetic ground states 
with incommensurate order, and more accurate methods 
to describe them in theory. 



Summarizing, we presented a successful formulation 
of the DMFT for strongly correlated magnetic systems, 
which opens a possibility of further applications in tran- 
sition metal oxides. In contrast to the earlier formula- 
tions based-Qn the modified second-order formula for the 
self-energy,E3 the present self-energy which describes the 
dynamical effects in the propagation of a hole coupled to 
spin fluctuations allows to obtain stable magnetic solu- 
tions: AF ordering at half-fiUing and SS in doped sys- 
tems. Although it is likely that. bett er variational states, 
possibly with stripe ordering,Q'E2rEil could be found it is 
expected that the presented spectral and optical proper- 
ties are generic for strongly-correlated systems with in- 
commensurate order parameter. A better understanding 
of the cuprates, however, requires a further development 
of theory which should be able to capture the gradual 
changes of local magnetic correlations in doped Mott- 
Hubbard systems under increasing temperature. 
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